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ABSTRACT 

Non-thermal electron populations are observed throughout the heliosphere. The relaxation of an electron 
beam is known to produce Langmuir waves which, in turn, may substantially modify the electron distribution 
function. As the Langmuir waves are refracted by background density gradients and as the solar and helio- 
spheric plasma density is naturally perturbed with various levels of inhomogeneity, the interaction of Langmuir 
waves with non-thermal electrons in inhomogeneous plasmas is an important topic. We investigate the role 
played by ambient density fluctuations on the beam-plasma relaxation, focusing on the effect of acceleration of 
beam electrons. The scattering of Langmuir waves off turbulent density fluctuations is modeled as a wavenum- 
ber diffusion process which is implemented in numerical simulations of the one-dimensional quasilinear kinetic 
equations describing the beam relaxation. The results show that a substantial number of beam electrons are ac- 
celerated when the diffusive time scale in wavenumber space To is of the order of the quasilinear time scale 
T q i, while when tq <k t cj i, the beam relaxation is suppressed. Plasma inhomogeneities are therefore an im- 
portant means of energy redistribution for waves and hence electrons, and so must be taken into account when 
interpreting, for example, hard X-ray or Type III emission from flare-accelerated electrons. 
Subject headings: Diffusion, Acceleration of particles, Sun:particle emission 
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1. INTRODUCTION 

The relaxation of an electron beam producing Langmuir 
waves is a fundamental process in space and laboratory plas- 
mas. As the Langmuir waves interact with the fast electrons, 
and their dynamics are strongly affected by variations in the 
ambient plasma density and as the heliospheric plasma den- 
sity is perturbed with different levels of inhomogeneity, the 
interaction of Langmuir waves with non-thermal electrons in 
turbulent inhomogeneous plasmas is an important topic in as- 
trophysics. 

In the solar atmosphere, the beam-plasma instability is 
known to be the cause of Type III radio bursts, with their 
characteristic rapid frequency drift being a result of the beam 
propagation through the decreasing density coronal plasma 
( Ginz burg & Zh elezniakovl 119581 : IZheleznvakov & Zaitsevl 
1970; Melrose 1990). Numerical simulations of the beam re- 
laxation in the inhomogeneous plasma of the solar wind can 
reproduce the broken power-law in the electron flux spectrum 
which i s observed near the E arth for impulsive solar electron 
events dReid & K ontar 2010). In solar flare environments, it 
has been shown that the beam relaxation is likely to cause 
spectral i ndex variations infe rred from hard X-ray observa- 
tions (e.g. Hannah et al. 2009). 

It was recognized already in the context of beam- 
plasma experiments (jBreizman & Rvutovl Il969t 

ILiperovskii & Tsvtovichl 119651) that me mere existence 
of an inhomogeneity in the ambient plasma density was 
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sufficient to produce a redistribution of Langmuir wave 
energy in wavenumber space. This idea was proposed in 
order to explain the common observations of electrons 
with energy above the injected one ( Breizma n & Rvutovl 
|1969|) in laboratory ex periments such as those carried out 
by Berezin et alj ( 11964ft on powerful relativistic beams. The 
interaction between the waves and the beam electrons is 
resonant, which implies that the spectral transfer of wave 
energy toward lower/higher wavenumbers allows interaction 
with faster/slower electrons respectively. Various mech- 
anisms may cause this spectral energy transfer, they may 
involve regular or random density variations, or non-linear 
interactions between Langmuir waves and low-frequency 
electrostatic compressive modes. For example, it is known 
that ion-sound waves can effectivel y scatter Langmuir waves 
to smaller wavenumbers (e.g., ILiperovskii & Tsytovichl 
[19651 iMelrosd [19741 lEscandd 119751: [PapadopouTosI 119751: 
lEscande & de Genouillac 1978). Low-freque ncy electromag- 
netic modes, such as ki netic Alfven waves (Bi an & Kontari 
120 10t iBian et all 12010) are compressive and thus their 
refractive power may also have a substantial effect on 
the spectral evolution of beam-driven Langmuir waves. 
Numerical simulations of the beam relaxation in a plasma 
with a constant d ensity gradien t were conducte d in (e.g. 
iKrasovskiil [l97l iKontarl l200lt iTsiklauril |2010|) and the 
combination of a regular density gradient with allowance for 
th e non-linear coupl i ng to i on sound waves was considered 
bvlKontar & Pecselil (120021) : I Yoon et alJ d2006h : lZiebell et alJ 
( 1201 11) . with the result that a very high level of Langmuir 
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wa ves was produ c ed at s mall wavenumbers. The simulations 
by iKontar et al.l (12012b found that a fluctuating density 
could lead to the appearance of accelerated electrons in a 
c ollisionally-rela x ing ele ctron beam. 

IVedenov et al.l (119671) first described the effect of ran- 
dom large scale density inhomogeneities as a diffusive 
transfer of Langmuir wave energy in wavenumber space. 
This d iffusion equation was used by Nishikaw a & Rvutovl 
d 19761) to describe the process of elastic scattering off ran- 
dom and time-independent density fluctuations, (see also 
iGoldman & Duboislll98l iMuschietti & Dumlll991l) . Elastic 
scattering results only in angular diffusion in wave-number 
space, and for the case of waves generated by beam electrons 
with a small angular spread, wave energy is transferred away 
from the region of excitation in £-space, which may lead to 
suppression of the beam-plasma instability. Neglecting the 
effect of wave reabsorption on the electrons, the role of the 
angular diffusion term has been accounted for in the kinetic 
equations describing the bea m-plasma system, by replacing 
it by a wave damping term (Muschietti et alj[1985t [Melrose 
1987). However, in general, energy conservation dictates that 
reabsorption of wave energy by the beam is accompanied by 
an acceleration of the electrons. In addition, inelastic scat- 
tering results in a change in the absolute value of the Lang- 
muir wavenumber, which may produce acceleration in the 
projected electron distribution. 

Here, we consider anew this important topic of beam re- 
laxation in a fluctuating plasma. In Section II, we derive a 
general expression for the diffusion coefficient in wavenum- 
ber space. In Section III, the wavenumber diffusion is im- 
plemented in numerical simulations of the one-dimensional 
kinetic equations describing the quasilinear evolution of the 
beam-plasma instability. We focus on the effect of acceler- 
ation of beam electrons in a fluctuating plasma and derive a 
condition for it to occur. In Section IV, we generalize the dis- 
cussion to 3D and discuss the role of both elastic and inelastic 
scattering in the acceleration of fast electrons. A summary of 
the results is presented in Section V. 

2. LANGMUIR WAVE SCATTERING AND DIFFUSION: ID CASE 

Let us start by considering a one-dimensional problem 
where the dynamics of both electrons and Langmuir waves 
are along x, the direction of the external magnetic field. The 
plasma density is written as n e [ \ +h(x, t)] with n e the constant 
background density and h(x, t) the relative density fluctuation, 
which is assumed to be weak, i.e. h(x, t) <k 1, a condition of- 
ten satisfied in the solar corona and the solar wind. 

In general, the characteristic wavenumber q of low- 
frequency density fluctuations is much smaller than the 
characteristic wavenumber k associated with high-frequency 
Langmuir waves, so we can make the WKB approximation 
and treat the Langmuir waves as quasi-particles. Ambient 
density fluctuations are associated with a change in the local 
refractive index experienced by the waves, and in the low fre- 
quency limit the Langmuir wave dynamics can thus be mod- 
eled b y Hamilton's equations of motion for t he quasi-particle s 
(e.g. lWhitham|[T96l IVedenov et all \l96% lZakharovl[l974l) . 
which are given by 



dk 
dt 
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dx 
dt 



where F{x,t) is the "refraction force" acting on the wave- 
packets, v g = 3Vj e k/a> pe is the group velocity and co pe is the 

local plasma frequency u pe = ^\nn e e 2 /m e . 

Equivalently, we can write a conservation relation describ- 
ing the evolution of the spectral energy density associated 
with Langmuir waves W(x, k, t) [erg crrT 2 ] 



dW(x,k,t) dW(x,k,t) , dW(x,k,t) 
• v„ F(x, ?)- 



dt 



dx 



dk 



0, (3) 



where W(x, k, t) is normalised to the energy density of Lang- 
muir waves E — J W(x, k, t)dk. 

From Equation ((TJ we see that random refraction induces 
stochastic change in the wavenumber, resulting in diffusion 
of the spectral energy density in £-space. An expression for 
the diffusion coefficient D is derived from the phase-space 
conservation equation Q following standa rd procedures (e.g. 
IVedenov & Velikhovll 9631: ISturrockll 1966b . as follows. 

The spectral energy density of Langmuir waves W(x, k, f) is 
decomposed into the sum of its average and fluctuating parts 
as W = (W) + W. Substituting this expression into Eq. © 
gives one equation for the average 

ot ok 
and one equation for the fluctuations 

dW dW d{W) 

where a term quadratic in the fluctuation amplitude has been 
neglected since the latter is assumed to be weak. Equation (0 
is integrated to give 

d(W)(x - v gT , k,t-r) 
dk ' 



W(x, k, t) 



f 

Jo 



drF(x - VgT, t — t) 



(6) 

which is substituted into Eq. © to give the equation describ- 
ing the diffusion of wave energy in £-space, 



d(W) _ d D d{W) 
dt ~ dk dk ' 



(7) 



where the diffusion coefficient D is expressed in terms of the 
auto-correlation function of the refraction force as 



D 



Jo 



dr(F(x, t)F(x - v g T, t - r)>. 



(8) 



The Fourier components F(q,£J) and the spectrum 
S f(q, Q.) of F(x, t) are defined through F (x, t) 
J dq J ^ d£lF(q, Q) exp [2m(qx — Qt)]dkd£l and 
(f(q,Q.)F(q',Cl')} = S F (q,Q.)6(q + q')6(Q. + £2'), re- 
spectively. As a consequence, the diffusion coefficient can be 
written in terms of S F (q, D.) as 

D=- I dq I dQS F (q,a)8(n-qVg). (9) 

Since the spectrum of the refraction force S F (q, Q) is re- 
lated to the spectrum of ambient plasma density fluctuations 
S n (q, Q) by S F (q, £2) = w pe (7rq) 2 S „(q, Q), the diffusion coef- 
ficient can finally be expressed as 

a> 2 n 2 r°° r°° 
(2) D=^f— dq\ dnq 2 S n (q,£l)5(£l-qv g ), (10) 

" t - kJ —CO *J —CO 
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where by definition (n ) = J dq J d£l S„(q, O). 

According to Eq. ( fTOl i. a quasi-particle with wave-number k 
interacts with a Fourier mode (Q, q) of the density fluctuation 
spectrum for which the resonance condition £2 = q\ g is satis- 
fied, which also means that (Q./q)(u> pe /k) - 3v 2 g . In the par- 
ticular case where density fluctuations are due to waves with 
a dispersion relation Q = Cl(q), then S„(q,D.) = S„(q)6(Q. - 
Q(q)) and D = (u 2 pe n 2 l2) jT dqq 2 S „(q)6(Q(q) ~ qv s ). 

For a Gaussian spectrum 



S„(q,Cl) = 



in 2 ) 
nqoQo 



exp 



9o 



(11) 



where qo and Qo are characteristic wavenumber and fre- 
quency, the diffusion coefficient reads 



D = ( J 3/2^0 2 
P V 



y 2\ 
V 



-3/2 



(12) 



where vo = Qo/qo- 

The diffusion coefficient depends on wave-number k 
through the group velocity v g = 3v 2 k/co pe . From Eq. (fT2l . 
we see that for a Gaussian spectrum there are two regimes of 
wavenumber diffusion. When the decorrelation velocity asso- 
ciated with ambient density fluctuations Vo is much larger than 
the group velocity of the Langmuir waves, i.e. when vo s> v g , 
then 



D 



(13) 



an expression which is independent of y g and hence k. Con- 
versely, when vo <s v g , then 

and the wavenumber diffusion becomes k dependent. 

It is also interesting to consider compressive fluctuations 
at a given characteristic frequency, so that S„(q,Q) = 
S„(q)6(Q. - Qo), with a power-law spectrum in wavenum- 
ber space wit h spectral index C > 1, as observed in the so- 
lar w ind (e.g. ICronyri[l97l iCelnikier et aDl!983b iRobinsonl 
1983]). Then 



S„(q,Q) = (n ) 



^|(-1U 9 



6(Q - Q ) (15) 



io j\qo) 

for q > qo and zero elsewhere. In this case we find that 

D = ^_1)(^) 3 V>. (16) 

Comparison of this last expression with Eq. ( fT4l i shows how 
the spectrum of density fluctuations affects the wavenumber 
dependence of the diffusion coefficient D. 

3. NUMERICAL RESULTS 

We now perform numerical simulations of the quasi- 
linear kinetic equations describing wave-particle inter- 
acti ons. These equati o ns are based on tho se give n 
bv IVedenov & Velikhovl d!963l) : IVedenov et ail d 19671) : 
iTsvtovichl dl995lT addin g a sponta n eous w ave emission 
term as in Zheleznvakov & Zaitsev (1970); Hannah et al. 



(2009 ): iReid et al.l (1201 ll) a nd collisional operators as in e.g 
iLifshitz & Pitaevskiildl98lT) : 
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(18) 

where F = 4ne 4 n e lnA/m 2 , with In A 20 the Coulomb 
logarithm. The last term in Eq. (U~8l describes the effect of 
wavenumber diffusion produced by the ambient density fluc- 
tuations, as discussed in the previous section. 

The initial electron distribution function /(v, f) [elec- 
trons crrT 3 (cm/s)" 1 ] is chosen as the superposition of a 
Maxwellian background and a Maxwellian beam: 



f(y,t = 0) = 



n e 



V2ttv 



■ exp 



Te 



,2 ^ 



«/, 



a/ttAv;, 



exp 



( (v-v,) 2 ) 



Av 2 



(19) 

and the initial spectral energy density of Langmuir waves 
W(k, t) [ergs cm" 2 ] is set to the thermal level: 



(20) 



Except where otherwise stated, the simulations below as- 
sume a background plasma with T e = = 1MK and 
/2n = 1GHz giving a background density of n e = 1.2 x 



3 . The beam parameters are nj, - 10 5 - 10 n e , 
x 10 9 cms _1 and Av;, = 0.3v/,. The spectrum of 



10 iu cm- 

V/, — 5 x luxms " anu zav^ 
background density fluctuations is taken first to be Gaussian, 
so the diffusion coefficient D is given by Eq. dT2l . We fix 
qo = 10~ 4 ko e . The magnitude of relative density fluctuations 

\/(« 2 ) ranges between 10~ 4 and 10" 1 , and Vo is taken between 
0.0 lv7> and Mj e . 

Figure (fl} is divided into four distinct panels which all show 
the evolution of the electron and the wave distribut ions as the 
intensity of the background density fluctuations V(" 2 ) is in- 
creased. In a homogeneous plasma, as shown in the top left 
panel, the electron beam is unstable to the generation of Lang- 
muir waves, which grow at a rate given by 



y ■ 



n b 



(21) 



The spectrum of Langmuir waves interacts with the beam 
electrons and causes them to diffuse in the resonant region in 
velocity space, given by (±> pe - kv, until the distribution func- 
tion reaches a plateau, on a time scale given by the quasilinear 
time: 

Tql = • (22) 

u) pe n b 

By time t = 100r f/ /, the beam has fully relaxed and the elec- 
tron distribution is flat between 6\Te (E = 1.5 keV) and 
16v7> (11 keV), as shown by the red lines in Figure Q] The 
collisional timescale, t co u = Vy IT and the timescale for 
collisional destruction of the plateau, t co //(v /v^) 3 ~ 1 s at 
v = 15v7> are both far longer than the quasilinear time. Col- 
lisional effects are therefore negligible. This remains true for 
the various values of beam density and background plasma 
parameters which will be considered below. 
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Fig. 1 . — Electron distribution /(v) (upper plots), and Langmuir wave spec- 
tral energy density W(k) (lower plots) for homogeneous (top left) and inho- 

mogeneous plasma, with qo = lCT 4 ^,,, vo/vjv = 0.3 and ^J(h 2 ) of 3.7xl0~ 4 
(top right), 3.7 X 1CT 3 (bottom left) and 1.2x 1CT 2 (bottom right). Lines from 
dark (blue in online version) to light show the beam relaxation during the first 
ten quasilinear times, while the thick line (red in online version) shows the 
state reached at r = IOOt,/. 

A finite level of background density fluctuations produces 
diffusion of the waves in fc-space which allows them to inter- 
act resonantly with electrons from a larger region of velocity 
space. In other words, the region of resonant wave-particle 
interaction is broadened due to the random wave refraction in- 
duced by the ambient density perturbations. This effect is well 
illustrated by the bottom left panel of Figure \T\ which shows 
the electron and wave dist ributions for a moderate level of 
density fluctuations, yj(n 2 ) = 3.7 x 1CT 3 . The Langmuir wave 
spectrum diffuses in wavenumber, widening the beam distri- 
bution and slowing the plateau formation. The electron dis- 
tribution /(v) is increased from the upper edge of the plateau 
to the highest velocities in the simulation, implying that elec- 
trons have been accelerated. In addition, the increased level of 
waves at large wavenumbers has a visible effect on the elec- 
tron distribution down to around 6vj> . 

The top right panel shows very weak density fluctuations, 

with y/(h 2 ) = 3.7 x 10~ 4 . The beam relaxes in a manner 
similar to the homogeneous case but the wave spectrum is 
broadened slightly in £-space, resulting in a wider plateau in 
the electron distribution. 

When the inhomogeneity becomes very strong, as illus- 
trated by the bottom right panel where yj(h 2 ) = 1.2 x 10~ 2 , 
diffusion still transports the waves out of their region of ex- 
citation in £-space but on a timescale which is much smaller 



Fig. 2. — The total energy of beam (left) and tail (middle) electrons, both 
normalised to the initial beam energy Eq, and the average energy of a tail 
electron (right), all evaluated at t = 100r 9 /, against the parameter R = td/t„i, 
which gives a measure of the degree of inhomogeneity. Beam and plasma 

parameters, and ranges for V(" 2 ) and vo are as given in the text. The asterisks 
mark the three levels of inhomogeneity shown in FigurefT] The dashed lines 
are polynomial fits to the data to illustrate its general trends. 

than their growth rate, so the wave level is barely increased 
above the thermal level. As we can see, by time t = lOOr,,/ 
the electron distribution remains essentially unchanged from 
the initial distribution. The main point is that for intense den- 
sity fluctuations, diffusive broadening of the Langmuir wave 
spectrum becomes large enough for the waves to be mainly 
reabsorbed by the thermal electrons at kAo e ~ L Since the 
energy density of waves generated is much smaller than the 
thermal energy of plasma, the effect on the distribution func- 
tion of thermal particles is negligible. This is the suppression 
regime of the beam-plasma instability which has often been 
discussed in the past, mainly in the context of the role of an- 
gular diffusion of wave energy induced by elastic scattering in 
3D. As shown here, the suppression of the beam-plasma insta- 
bility can also occur as a result of inelastic scattering. Here, 
we focus on the effect of electron acceleration due to wave re- 
absorption by the beam electrons as opposed to reabsorption 
by the the thermal electrons. A discussion of electron acceler- 
ation in 3D, as a result of both elastic and inelastic scattering, 
is presented in the next section. 

Clearly, the main effect of increasing the level of density 
fluctuations is to increase the rate at which wave energy is 
transported in fc-space. The time scale associated with such 
diffusive process is 

(Ak b ) 2 
D(k b ) 



where k b = a) pe /Vb an d A^£> = oi pe (AVi,/v?) , with parameters 
to pe ,Vb and Av/, as given above. Therefore, it is natural to 
define a non-dimensional number R, which is the ratio of the 
diffusion time t/j and the quasilinear time r 9 ;, i.e. 



Tql 



(24) 



To quantify things further, let us define a "beam" region in 
velocity space that includes all electrons with velocity above 
6vt>, and a "tail" region above 16vt>. The energy of the ini- 
tial Maxwellian beam is given by Eq. We can then deter- 
mine the total energy in the beam and tail electrons at time 
t = 100r fy /, which indicates the extent to which the beam re- 
laxation is suppressed and electrons are accelerated. In the 
homogeneous case, E\, eam = Q.I6E0 and E ta u = Q.2E^ respec- 
tively. In Figure |2l we plot the ratios Ei, eam /Eo, E tai ilE[) and 
Etaii/ntaii at t = 100t 9 / as functions of the parameter R, for a 

range of values of V(" 2 ) and vo. The three cases in Figure Q] 
are marked by asterisks, and represent three distinct regions in 

R. In the strong inhomogeneity case, with ^{n 2 ) = 1.2 X 1CT 3 
and R <sc 1, the beam-plasma instability is suppressed, so 
the beam remains close to its initial Maxwellian form and 
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Fig. 3. — As Fig[2] the total energy of beam (left) and tail (middle) electrons, 
both normalised to the initial beam energy Eq, and the average energy of a 
tail electron (right), all alt = lOOr,; against R = Tohql- Simulation results 
for various values of the beam, plasma and density fluctuation parameters 
are shown. Ranges for V("), vo are as m me text - For 0J pe = 1GHz we 
show: Gaussian fluctuations with m, = 10 5 (crosses) and nt = 10 4 (squares), 
powerlaw fluctuations with indices | and 5 (triangles). Finally for cj pe = 
200MHz we show Gaussian fluctuations with nt = 10 3 (diamonds). 

we find Ebeam = 0.99Eo and E ta u = 0.21.Eo. Very weak 
inhomogeneity, with ^j{h 2 ) = 3.7 x 10~ 4 corresponding to 
R » 1, has little effect on the beam relaxation, so the beam 
energy is close to the homogeneous value, but the tail en- 
ergy is slightly increased due to the broadened plateau, giving 
Ebeam = OJ8E0 and E la \\ = Q.23Eo The intermediate case, 
with V^ 5 ) = 3.7 x 1CT 3 and R ~ 1, gives E beam = 0.95E 
and E ta ii = 0.35£o- Both the total tail electron energy (middle 
panel) and the average energy of a tail electron (right panel), 
show a significant increase over more than two orders of mag- 
nitude in R around R ~ 1. To summarize we have identi- 
fied three distinct regimes of the beam plasma instability in a 
fluctuating plasma which are controlled by the parameter R. 
When R » 1, the density fluctuations have a weak refrac- 
tive power and the relaxation proceeds as in a homogeneous 
plasma. When R <k 1, diffusive broadening of the wave spec- 
trum is large, the waves are mainly reabsorbed by the ther- 
mal particles, their level stays small and the instability is sup- 
pressed. When R ~ 1 diffusive broadening of the wave spec- 
trum is such that a substantial part of the wave energy is now 
also reabsorbed by the beam electrons leading to acceleration 
of the beam particles. Using Eq. (1121 1. we may express the 
condition for strong acceleration, R ~ 1 as 



n b I Av, 



2 2 



V 2 ^ 
V 0J 



-3/2 



(25) 



As discussed in Section 2, for a Gaussian spectrum of den- 
sity fluctuations there are two extreme regimes of wavenum- 
ber diffusion, where the characteristic group velocity of the 
beam driven Langmuir waves, v s is much larger or smaller 
than Vq. Let us notice that all of the values for Vq used in the 
above simulations lie in the transitional region v g ~ vo. The 
effect of vo on the ^-dependence of the diffusion coefficient 
and hence on tq is therefore also important, and accounts for 
the vertical scatter in the points in Figure|2] 

Figure [3] shows the total beam and tail energies and av- 
erage tail electron energy for several values of background 
density and beam density, for both Gaussian and power-law 
density fluctuation spectra. The distribution of the beam and 
tail total energies with respect to the parameter R is essen- 
tially unchanged from Figure [2] We note that collisional ef- 
fects remain unimportant for all parameters chosen. Lower 
beam densities than those shown produce insufficient levels 
of Langmuir waves to cause acceleration, and the parameter 
R is unaffected by changes in local plasma frequency. The 
results for power-law fluctuation spectra are similar to those 
for a Gaussian spectrum. We may therefore expect significant 



electron acceleration to occur whenever R ~ 1, regardless of 
the specific spectrum of fluctuations, or the specific values of 
the other parameters. 

4. GENERALIZATION TO 3D 

Using the same methodology developed for the one- 
dimensional case, we now generalize the diffusive descrip- 
tion of Langmuir waves to three-dimensions. The main dif- 
ference with the ID case is that in 3D refraction can change 
both the absolute value and the orientation of the wavevector 
k. The conservation equation for the spectral energy density 

W(x, k, t) [erg], is now 



dW(x, k, t) 
di 



+ v s .Vff(i, t, t) - -w pe Vh(x, t) 



dW(x,k,t) 
8k 



0. 

(26) 



where yj, = (3\j e /a) pe )k is the group velocity. In 3D the dif- 
fusion equation takes the form 



dW(k, t) 



d_ D dW(k, t) 
dk~i U dkj 



(27) 



with the diffusion tensor given by 

-/3, 



D u (k) = 2naj 2 pe J ' dD. J q^jSM ^ {^-q- v* g ) , 

(28) 

where S„{q, Q) is the spectrum of the density fluctuations and 
J c/kW(k, t) = E. In the particular case where density fluctu- 
ations are due to waves with a dispersion relation Q = Cl(q), 
then S„(q, Q) = S„(q)6(D. - D.(q)) and 

D u (k) = 2noj 2 pe Jd^q q iqj S n (q)8 - 4 ■ v^) . (29) 

In spherical coordinates with k = (k, 9, <p), the general equa- 
tion for diffusion (assuming azimuthal symmetry) is 



dW(k,t) _ I 1 d I 2 d_\ J 



S -(Sb))"** 

(30) 

with coefficients given in the Appendix below by Equation 
Ml 

Expressions for the only non-zero components of the diffu- 
sion tensor Da and Dgg are also derived in the Appendix lAl 
and are given by 



D, 



0)„ 



2l67Tvl e (kA de ) 



89 



2An\ Te kA de 



J dqq 3 1 - ^ 



dqn(qfqS n (q) (31) 



2\ 



1-1 m ) 

3v T ekA de q) t 



S„(q) (32) 



Note here that due to the resonance condition Q(g) < 
3vrekA de q the second expression is always positive. 

The particular case of elastic scattering is recovered by set- 
ting Q. = in Equations OTb and d32i l in which case we find 
that 

D kk = (33) 



and 



2\it\j e kA, 



de 



f 



dqq S n (q), 



(34) 
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as given by Mu schietti et alj (|T985). Diffusion then occurs in 
angle only at a rate independent of 6 and tends to isotropise 
the Langmuir wave spectrum. 

It is instructive to rewrite the diffusion equation in terms of 
the independent variables k and k\\, where k\\ is the component 
of the wavevector parallel to the beam direction, given by k\\ = 
kcosO. In the limit cos# ~ 1, or equivalently k ~ k\\, we 
obtain 

| W( *"< ?) = |^ W( ^ " 2 ^ W(kh tY 

(35) 

This shows that both magnitude and angular change in wave- 
vector, as given by D kk and Dgg, contribute to the redistribu- 
tion of wave energy in wavenumber parage/ to the beam. This 
transport of wave energy is diffusive for the former and con- 
vective for the latter. In both cases, waves shifted to smaller 
parallel wavenumbers can be absorbed by beam electrons at 
larger parallel velocities resulting in acceleration of tail elec- 
trons. A large body o f work has been devoted to t he role of an- 
gular diffusion (e.g. Nishikaw a & Ryutovll 19761) with respect 
to the suppression of the beam-plasma instability, mainly in 
the context of elastic scattering, i.e. D kk = 0. The electron 
distribution was assumed to be fixed and the effect of scatter- 
ing of waves by density fluctuations was modeled solely as 
wave damping. However, the energy lost by the waves due 
to damping is transferred to the electrons leading to electron 
acceleration. 

A full numerical treatment of this 3D diffusion and its accel- 
eration effects is not possible using our ID quasilinear code, 
and projecting the 3D case into ID is uninformative. How- 
ever, as seen in the previous section, acceleration of beam 
electrons in the ID case occurs when the transport time scale 
in £-space is of the order of the the quasilinear time scale, 
for a wide range of beam, plasma and fluctuation parame- 
ters. Angular diffusion convects wave energy in k\\, always 
toward small h and on a time scale given by Tgg - k^/Dgg 
(see Equation i35[). Thus we expect angular diffusion to also 
lead to electron acceleration, confirmed by, for example, the 



PIC simulations of iKarlickv & Kontarl d2012). which include 
the effect of long wavelength density fluctuations on Lang- 
muir waves, and find such an acceleration effect. From our 
diffusion treatment we may estimate that this effect will be 
significant when Tgg ~ T q i, whether scattering is elastic or not. 
It is also natural to expect that a necessary condition for the 
electron distribution function to remain fixed, and hence the 
beam-plasma relaxation to be suppressed, is that Tgg <K T q i 
otherwise when Tgg ~ T q i the effect of electron acceleration 
studied in this work needs also to be considered in the overall 
wave-particle energy budget. 

5. SUMMARY 

In summary, we have considered the effects of the diffusion 
of Langmuir waves in wavenumber space numerically, using 
a 1 -dimensional model. We found the potential to suppress 
the beam-plasma instability when the diffusion is sufficiently 
fast, but also the possibility of a significant acceleration of 
beam electrons. The transition from acceleration to suppres- 
sion is controlled by the ratio of the quasilinear and £-space 
diffusive timescales, with the most efficient acceleration oc- 
curring when the quasilinear time is close to the timescale 
for diffusion of Langmuir waves in wavenumber space. This 
relationship was found to hold for a wide range of beam, 
plasma, and fluctuation parameters. In addition, while not 
presented here, previous work found a similar acceleration 
effect to o ccur during the co llisional relaxation of a power 
law beam (Kon tar et alj|20 12). We have also derived the dif- 
fusion coefficients for the general case of inelastic scattering 
in 3-dimensions. This constitutes an extension of the elastic 
scattering theory previously discussed in the literature. A co- 
efficient D kk " and a correction to the coefficient Dgg pro- 
portional to Q. 2 were found. From the diffusion equation and 
the effects of diffusion in ID, we infer that scattering in 3D 
will be capable of producing electron acceleration, provided 
the transport time-scale in wave-vector space is of the order of 
the quasilinear timescale, while the beam plasma relaxation is 
suppressed only under the condition that the former is much 
smaller than the later. 



APPENDIX 

CARTESIAN AND SPHERICAL REPRESENTATIONS OF THE DIFFUSION TENSOR 

The diffusion equation in Cartesian coordinates is 

d,W(t t) = d ki Dijd kj W{t t) (Al) 

with coefficient given by 

D u (k) = 2noj% fdClj qiqjSM Q)6 (Q - 4 • vj) (A2) 
where S„(q, Q) is the spectrum of the density fluctuations. 

Definition of coordinates 

In Cartesian coordinates, we define one axis to be parallel to the beam direction, labelled as ||, and have two mutually perpen- 
dicular axes labelled, ±\, ±2, giving a standard right-handed Cartesian coordinate system. We also define a spherical coordinate 
system, with 6 the angle to the beam direction, and <p the azimuth, measured clockwise around the beam direction. 

In these spherical coordinates we write the Langmuir wavevector as k = (k, 8, <p). Assuming azimuthal symmetry, we may 

define the coordinates such that the azimuth of k is zero. The diffusion tensor can then be transformed from the Cartesian 
coordinates, in which it has components D,j (D\\\\, D\\ ±1 , D||x 2 etc) as given by Equation lA2l into these spherical coordinates. The 
diffusion equation may also be transformed, becoming 

d ki D u d kj W(k, t) = Updk (k 2 D kk d k + kD ke dg) + -^dg smO^dg + ^yd k H W& t) (A3) 
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There is no diffusion in <p due to azimuthal symmetry. 
By analogy with the Cartesian expression in Equation lA2l we define the following quantities 

qm = 1± 1 <1±, cos2 — 2<3 , j_ 1 gu sin 6 cos 6 + q^ sin 2 6 (A4) 
qkk = q±i q±i sin 2 6 + 2q ±i q\\ sin 6 cos 6 + ^h^h cos 2 9 (A5) 

and 

qw = qek = sin cos 0(q ±l q ±1 - q\\q\\) + (cos 2 6 - sin 2 0)q M q\\ (A6) 

which allow us to write simply 

Dijit) = 2no? pe j dQ, J J^dq j dp J d0 q u S (fl - 4 ■ (A7) 

where j - 9, k, (p. 

Expressions for the diffusion coefficients 

Assum ing a dispersion relation Q = Q(q) (the general case being intractable) we write S„{q,Q) = S„(q)S(Q. - Cl(q)) so that 
Equation l A7 I becomes 

Dy$ = 2ncj 2 pe J J^dq jdpjdj> qijS(q)d(n -q-v g ) (A8) 

Substituting 

k-q = kq({\ -yu 2 ) 1/2 (l -p 2 ) 1/2 cos 4> + (lp) (A9) 

into the delta function, and integrating over (f>, the resonance condition becomes 

£1' — p.p 



(\-^H\-p 2 ) 



1/2 



(A10) 



where £l'(q) = Q.(q)u> pe /(3v 2 e kq) (resona nce a lso requires O' < 1). 
We use this to evaluate Equations I A4I to |A6l finding 

qw = 



q kk = q 2 n'(q) 2 
: q 2 [p 2 + fl 2 Q.'(q) 2 - 2tfW( 9 )) /(l - fl 2 ) 



The diffusion coefficients are then 

,3 



"Te 

and 

,3 



D kk = - J dqq 3 J dfxQ.'(qfS n (q) ((1 - /i 2 )(l - p 2 ) - O' 2 + 2^Q' - fi 2 p 2 )~' /2 (Al 1) 



Uaa — 



^ J dqq* j dp * + ^ ^-^% ((lVXl-^-fl^Wy-^r (M2) 

where the limits on the p integrals are the solutions of (1 -/U 2 )(l —p 2 )-Q. ' 2 +2p.p£2' -p 2 p 2 = i.e /i ± = yuQ' + (l -// 2 )^ 2 (1 - Q' 2 ) 1 ^ 2 
(noting that O' < 1). 

Assuming isotropic fluctuations, the p integrals can be calculated, and the results are 

°kk = — f dqqn(q) 2 S„(q) (A13) 



2l6nvl ikXde) 



S»(9). (A14) 



The special case of elastic scattering, that was treated by Nishikaw a & Rvutovl dl976l) : lMuschietti et al.l (fl985) may be simply 
recovered by taking £2 = 0, for which we find 

D kk = (A15) 
D ee = j dqq 3 S„{q) (A16) 
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